% Written by: Orhan Torul
% Istanbul, Turkey, May 2017
function Gvec = sgushimer(varshimer,fy1,fy,fx1,fx)

global ga ha;

a11_g=varshimer(1);
a12_g=varshimer(2);
a21_g=varshimer(3);
a22_g=varshimer(4);
a31_g=varshimer(5);
a32_g=varshimer(6);
a41_g=varshimer(7);
a42_g=varshimer(8);
a51_g=varshimer(9);
a52_g=varshimer(10);
a61_g=varshimer(11);
a62_g=varshimer(12);
a71_g=varshimer(13);
a72_g=varshimer(14);
a81_g=varshimer(15);
a82_g=varshimer(16);
a11_h=varshimer(17);
a12_h=varshimer(18);
a21_h=varshimer(19);
a22_h=varshimer(20);
f_y_1=fy1;
f_y=fy;
f_x_1=fx1;
f_x=fx;
ga=[a11_g a12_g; a21_g a22_g; a31_g a32_g; a41_g a42_g; a51_g a52_g; a61_g a62_g; a71_g a72_g; a81_g a82_g];
ha=[a11_h a12_h; a21_h a22_h];

Gara =f_y_1*ga*ha + f_y*ga + f_x_1*ha + f_x;
Gvec=Gara(:);
